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Systems living in complex non equilibrated environments often exhibit subdiffusion characterized 
by a sublinear power-law scaling of the mean square displacement. One of the most common models 
to describe such subdiffusive dynamics is the continuous time random walk (CTRW). Stochastic 
trajectories of a CTRW can be described in terms of the subordination of a normal diffusive process 
by an inverse Levy-stable process. Here, we propose an equivalent Langevin formulation of a force- 
free CTRW without subordination. By introducing a new type of non-Gaussian noise, we are able to 
express the CTRW dynamics in terms of a single Langevin equation in physical time with additive 
noise. We derive the full multi-point statistics of this noise and compare it with the scaled Brownian 
motion (SBM), an alternative stochastic model describing subdiffusive dynamics. Interestingly, 
these two noises are identical up to the 2nd order correlation functions, but different in the higher 
order statistics. We extend our formalism to general waiting times distributions and force fields 
and compare our results with those of SBM. In the presence of external forces, our proposed noise 
generates a new class of stochastic processes, resembling a CTRW but with forces acting at all times. 


I. INTRODUCTION 

Many systems in nature live in complex non- 
equilibrated or highly crowded environments, thus ex¬ 
hibiting anomalous diffusive patterns, which deviate from 
the well known Fick’s law of purely thermalized systems 
w- Their distinctive feature is the power-law scaling 
of the mean-square displacement (MSD) PJj3]: 

((Y(t)-Y 0 ) 2 )~ t «, (l) 

where (•) indicates the ensemble average over different 
realizations of the stochastic process Y (t) describing the 
dynamics, usually either a velocity or a position, Y 0 is 
its initial condition and a £ R + . While Fick’s law is 
recovered by setting a = 1, thus predicting for normal 
diffusion the typical liner scaling of the MSD, we can dis¬ 
tinguish between different types of anomalous behaviour. 
Indeed, we define subdiffusion if 0 < a < 1 and superdif¬ 
fusion if a > 1, which correspond to processes dispersing 
with a slower or faster pace than Brownian motion re¬ 
spectively. 

Examples of such anomalous processes were first found 
in physical systems, such as charge carriers moving 
in amorphous semiconductors, particles being trans¬ 
ported on fractal geometries or diffusing in turbulent flu¬ 
ids/plasma or in heterogeneous rocks (see [2] and refer¬ 
ences therein). However, with the recent improvements 
of experimental techniques in biology, joint position- 
velocity datasets have been obtained, which are revealing 
the existence of many more examples in living systems. 
On the one hand, cell migration experiments have re¬ 
vealed a characteristic superdiffusive scaling of the MSD 
along with many more features deviating from standard 


* Correspondence to: a.baule@qmul.ac.uk 


Brownian models, e.g., non Gaussian probability density 
functions (PDFs) for the position and/or the velocity of 
the moving cell and power-law long time scaling of the 
velocity auto-correlation functions (see Refs. [5] and ref¬ 
erences therein). On the other hand, biological macro- 
molecules and/or organelles often exhibit a subdiffusive 
scaling of the MSD, while moving in molecular crowded 
environments. These can be prepared ad hoc for in vitro 
experiments, e.g., by using solutions of surfactant mi¬ 
celles [7] or polymer networks [8] to name just a few tech¬ 
niques, or found in vivo, e.g., in the cytoplasm [9] or the 
cells’ membrane uni, whose viscoelastic properties have 
recently been found to play a major role in determining 
the anomalous diffusion mm- Furthermore, biological 
systems may also exhibit a rich dynamical behavior, such 
as non linear MSDs, showing crossover between different 
scaling regimes |T3Hl9] at different timescales, and/or a 
dependence of the corresponding diffusion coefficients on 
energy-driven active mechanisms mum. We refer 
the interested reader to |21| for a recent review on anoma¬ 
lous diffusive systems. 

Considering this wide, though not exhaustive, variety 
of different anomalous behaviours, one needs to have a 
tool-kit of well studied models with which trying to fit 
the experimental data and infer the specific microscopic 
processes underlying the observed dynamics. Here we 
focus on subdiffusive processes, for which many models 
have been introduced so far, which are capable of repro¬ 
ducing the characteristic scaling of Eq. ([!]), while still 
showing distinct features if we look at other properties, 
like the multipoint correlation functions |22H25] . Among 
the most commonly applied to data analysis, we find the 
continuous-time random walk (CTRW) |2J [25] and the 
scaled Brownian motion (SBM) [271 - I3Q] . 

In the seminal paper [25], the CTRW was introduced 
as a natural generalization of a random walk on a lat¬ 
tice, with waiting times between the jumps and their 
size being sampled from general and independent prob- 
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ability distributions. Only later, a convenient stochastic 
representation of these processes was derived in terms 
of subordinated Langevin equations EH> which provided 
a suitable formalism to derive their multipoint correla¬ 
tion functions :22] (33]. Although the focus was first on 
power-law distributed waiting times, which indeed pro¬ 
vided Eq. 0 exactly for all times, recent works adopted 
more general distributions [32H33] . thus being able to 
model the crossover phenomena so often occurring in bi¬ 
ological experiments. 


On the other hand, the SBM has been recently intro¬ 
duced as a Gaussian model of anomalous dynamics [27] . 
providing the same scaling of Eq. 0 for all its dynam¬ 
ical evolution. If B(t) is a usual BM, its scaled version 
is defined by making a power-law change of time with 
exponent a: B(t a ). Although being commonly used to 
fit data [7, |2H [35], it has recently been shown to be a 
non stationary process with paradoxical behaviour under 
confinement, i.e., in the presence of a linear viscous-like 
force, as the corresponding MSD unboundedly decreases 
towards zero. This is suggested to be ultimately caused 
by the time dependence of the environment, either of the 
temperature or of the viscosity. As a consequence, it has 
been ruled out as a possible alternative model of anoma¬ 
lous thermalized processes [5Pj . 


In this paper, we derive a new type of noise, which 
allows us to express a free diffusive CTRW in terms of 
a single Langevin equation in physical time. We pro¬ 
vide the full characterization of its multipoint correlation 
functions and we compare them with those of the noise 
driving a SBM. Both purely power-law waiting times and 
general waiting time distributions are discussed [331 EH| ■ 
We find that the correlation functions are identical up to 
the two point ones, but different for higher orders: the 
noise driving SBM is a Gaussian noise, while our new 
noise driving a CTRW is clearly non-Gaussian. Here, all 
odd correlation functions vanish, as for Gaussian noise, 
but the even ones do not satisfy Wick’s theorem [5B]. The 
newly defined noise enables us to define a class of CTRW 
like processes with forces acting for all times, which are 
different from the corresponding standard CTRWs. Fur¬ 
thermore, we revisit the behaviour of the SBM under con¬ 
finement and show that its MSD correctly converges to a 
plateau as it is typical of confined motion EZ], provided 
that we use more general time changes with truncated 
power-law tails. This suggest that the anomaly observed 
in [3D] is mainly due to the localizing effect of the exter¬ 
nal linear force, which is able to trap the particle in the 
zero position if we allow for infinitely long waiting times 
between the jumps to eventually occur in the long time 
limit. 


II. CTRWS, SCALED BM AND 
GENERALIZATION TO ARBITRARY WAITING 
TIMES’ DISTRIBUTIONS AND TIME 
TRANSFORMATIONS 

We recall in this first section definitions and proper¬ 
ties of the free diffusive CTRW and SBM, which will be 
useful later in the discussion. We are mainly interested 
in their stochastic Langevin formulation and in both 
their Fokker-Planck (FP) equation and MSD. We then 
generalize these results to the case of arbitrary waiting 
times’ distributions or time transformations for CTRW 
and SBM respectively. 

Throughout the discussion, / indicates the Laplace 
transform of a function f(t) defined on the positive half 
line: /(A) = C {/(£)} (A) = / 0 + °° e~ xt f(t) dt. Moreover, 
ipi * ifi 2 denotes the convolution of two functions f\ and 
ip 2 , defined on the positive half line: (ip\ * f 2 )(t) = 
f 0 fi(t — t)<p 2 (t) dr. We remark that the correspond¬ 
ing definitions for functions of multiple variables follow 
straightforwardly. 

A. CTRW 

A Langevin representation of a CTRW was first pro¬ 
posed in 31], where the method of the stochastic time- 
change of a continuous-time process is used. Its set-up 
consists in introducing two auxiliary processes A'(s) and 
T(s), which we assume for now to be purely diffusive and 
Levy stable with parameter a (0 < a < 1) respectively. 
They both depend on the arbitrary continuous parame¬ 
ter s and have dynamics described in terms of Langevin 
equations [311 : 

X(s) = V2a£(s) (2a) 

T(s) = V(s) (2b) 

where £(s) and r](s) are two independent noises. For A'(s) 
to be a normal diffusion, we require £(s) to be a white 
Gaussian noise with (£(s)) = 0 and (£(si)£(s 2 )) = S(s 2 — 
s i). On the other hand, r](s) is a stable Levy noise with 
parameter a (0 < a < 1) [33]. The anomalous CTRW is 
then derived by making a randomization of time, i.e., by 
considering the time-changed (or subordinated) process: 
Y(t) = X(S(t)), with S(t) being the inverse of T(s), 
defined as a collection of first passage times: 

S(t) = mf{s:T(s)>t}. (3) 

s>0 

The process Y ( t ) is easily shown to satisfy Eq. ([!]) exactly 
for all its time evolution, by recalling that the probability 
density function (PDF) of S(t) has the Laplace transform 
h(s,X) = A“ -1 e -sAQ [22J and that (A' 2 (s)) = 2crs. In¬ 
deed, we obtain in Laplace space: 

(y 2 (A))= J + °° ( X \s))h(s,X)ds=^, (4) 
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whose inverse transform confirms its anomalous scaling: 

< y2 <*» = WT^f- < 5 > 

As expected, this same MSD is obtained by taking the 
diffusive limit of the microscopic random walk formula¬ 
tion of the CTRW, where we allow for asymptotically 
power-law distributed waiting times between the jumps 
of the walker, whose sizes are drawn from a distribution 
with finite variance [2|. In this limit, the model also pro¬ 
vides a fractional diffusion equation for the PDF of F( t ): 




( 6 ) 


where D a is a generalized diffusion coefficient and 
V\- a m = rp|/o(i-' r )“" 1 /W is the Riemann- 
Liouville time-integral operator, which makes the non 
Markovian character of the CTRW evident. It is then 
natural to investigate if the set of Eqs. ( 2apb ) can 
give this same FP equation. This has been proved in 
[34, 39} HD], with the specification: D a = T ^i +a ) , thus 
confirming the equivalence in the diffusive limit of the 
original random walk model and of the subordinated 
Langevin Eqs. ( 2apb ). 

We remark that the formulation of CTRWs as a sub¬ 
ordinated normal diffusive processes can be considered 
as the continuum limit of the original renewal picture of 
Montroll and Weiss [28j. Here, the position Y(t) of a 
CTRW is characterized by two sets of random variables 
{(&; r li)}i,...,N(t)> with N(t) being the number of jumps 
up to the time t. The random variables £» and rji spec¬ 
ify the amplitude of the jumps occurring at the random 
time tj, i.e., £* = F(tj) — Y(ti- 1 ), and the waiting times 

between two successive jumps, i.e., rji = U — _i. Thus, 

Y (t) is obtained by summing all the variables 


N(t) 

Y{t) = Y,Zi- ( 7 ) 

*=l 

In the present discussion, we assume {^} and { 77 ^} sep¬ 
arately to be i.i.cl random variables and each to be 
independent of rji. On the other hand, we obtain by di¬ 
rect integration of Eqs. pipbfr 

Y(t) = X(S(t)) 

X{t) dr 

... fS(t) 

= y/2a / £(t) dr. (8) 

Jo 

Therefore, Fogedby’s approach m describes the result¬ 
ing trajectory of the random walk in the continuum limit 
by parametrizing both the path of the walker X(-) and 
the time elapsed T(-) with an arbitrary continuous arc- 
length s. The stochastic process S(-) is the inverse of T(-) 



and measures the arc-length as a function of the physi¬ 
cal time. S(t) thus represents the continuum limit of the 
random variable N(t) that counts the number of steps in 
the renewal picture. 


B. SCALED BM 

If instead of a stochastic time change, we consider the 
deterministic time transformation t —► t* = t a in the 
normal diffusive process X (t) (now in the physical time 
t), we obtain the SBM: F*(t) = X(t*). Its equivalent 
Langevin equation is given by [27HS3: 

Y*(t) = V2aat a ~ 1 £(f), (9) 

with £(t) being a white Gaussian noise (with the same 
properties as before, but in the physical time t). By using 
Eq. ([9]) we can prove straightforwardly that the MSD of 
F* (f) is the same as Eq. © and that the corresponding 
FP equation is given by: 

-^P(y, t) = a <7 l Q y - 2 r (y^ *)> ( 10 ) 

which has time dependent diffusion coefficient [41]. This 
process preserves all the properties of Brownian motion 
it is indeed Gaussian with time-dependent variance 
and Markov, as the monotonicity of the time change pre¬ 
serves the ordering of time. Furthermore, F*(t) is self¬ 
similar and it has independent increments for non over¬ 
lapping intervals. However, differently from Brownian 
motion, it is strongly non stationary [30l . Furthermore, 
F* (t) turns out to be the mean-field approximation of the 
CTRW, as it describes the motion of a cloud of random 
walkers performing CTRW motion in the limit of a large 
number of walkers 129] . Recent investigation have also 
shown that SBM exhibits rich aging properties, which 
strongly differentiates it from the standard BM [12] . 


C. ARBITRARY WAITING TIMES’ 
DISTRIBUTION AND TIME 
TRANSFORMATIONS 


In this section, we first focus on the generalization of 
Eqs. ( 2a||2b ) to arbitrary waiting time distributions of 
the underlying random walk [32ll34l l43l . This extension 
is obtained naturally by choosing a different process T(s) 
with the only assumption of it being non decreasing in 
order to preserve the causality of time. Thus, we consider 
r)(s) in Eq. (2b I to be an increasing Levy noise with paths 
of finite variation and characteristic functional 1441: 


G[k(r)] = (e~ fc O ?) ( T ) dT ^ = e -/ 0 +DO $ ( fe ( T )) dT . 

( 11 ) 

Here $(fc(r)) is a non negative function with $(0) = 0 
and strictly monotone first derivative, while k(r) is a test 
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function. We recall that for 4>(s) = s“ we recover the 
CTRW model. Under these assumptions, the integrated 
process T(s) is a a one-sided increasing Levy process with 
finite variation. Furthermore, we assume r](s) to be in¬ 
dependent on the realizations of £(s) in Eq. |2a| . As a 
consequence of the finite variation and the monotonicity 
of the paths of T(s) respectively, S(t) has continuous and 
monotone paths, with this second property implying the 
fundamental relation |22j : 


= l-e(i-T(s)). (12) 


Similarly to Eq. Q, we can derive the corresponding 
MSD by recalling that h(s, A) = 5^d e - s$ ( A ) [33, 34]: 


{Y\t))=2a [* K{t) dr, 
Jo 


(13) 


for K{t) being related in Laplace space to <f>(s) by: 

i<W = jE- (14) 



Furthermore, the PDF of Y (t) is obtained by solving the 
generalized FP equation [34! : 

§t P{v ' t) = , ’W§il K(t-OP(v.r)dr, (15) 

whose solution in this particular case can be found for 
general 4>(s) in Laplace space: 


P(y, A) = 



(16) 


We look as an example at the case of a tempered stable 
Levy noise with tempering index /x and stability index a 
[45]. which is obtained by setting 4>(A) = (/z + A)“ — 
/z“, i.e., K(t) = fa-1 E ata ((nt) a ) [36]. As already 
pointed out, the CTRW case is recovered by setting /i = 
0, meaning that we do not truncate the long tails of the 
distribution, thus accounting for very long waiting times 
with a power-law decaying probability of occurrence. We 
plot in Figure[l]the numerical Laplace inverse of Eq. (16) 
(main) and the corresponding MSD (inset) at a fixed time 
t = 1000 (dotted line in the inset), which is given by [47] : 


(Y 2 (t)) 


2 a_ _ 1 + y^ 7 (/A; an) 
u a r (an) 

r L n—0 v ' 


(17) 


with "f(x; a) = Jq e t t a 1 being an incomplete gamma 
function, leading to the asymptotic behaviour |34j [47] : 


Figure 1. (Colors online) PDF (main) and MSD normalized 
to t (inset) of an anomalous process Y(t) obtained by sub¬ 
ordination of a pure diffusive process with a tempered stable 
Levy process of tempering index fx and stability parameter 
a = 0.2. The PDF is obtained by numerical Laplace inversion 
of Eq. ( |16| ) at t = 1000 (black dotted lines in the inset) [48j. 
The smooth transition from the non-Gaussian PDF typical 
of CTRWs (fi = 0) and the Gaussian one of normal diffusion 
(fi —» Too) is evident, along with the corresponding transition 
from anomalous to normal scaling of the MSD for increasing 
fi at a fixed time. Simulations, obtained with the algorithms 
of l49l l50j . agree perfectly with the analytical results. 


subdiffusive to normal behaviour, but their MSD scales 
as a power-law for all times. As expected, for /i = 0 we 
recover the typical non Gaussian shape of the PDF of a 
free diffusive CTRW [2]. However, for increasing values 
of fi, the PDF of Y(t), although still being non Gaussian, 
broadens, thus getting closer to a Gaussian. This has also 
evident consequences on the dynamical behaviour of the 
MSD, which for increasing values of fi goes from a pure 
subdiffusive scaling to a normal one (inset). 

We now discuss the corresponding extension of the 
SBM to arbitrary time transformations involving the 
kernel K(t ) obtained by Laplace inverse transform of 
Eq. ( |l4[ . We then generalize Eq. by adopting K(t) 
as the time dependent coefficient of the white noise: 

Y{t) = ^K(t)m = at), (is) 


<r 2 (t)> 


2 <7 

r(l+a) 

'v <y ' 


t a 
“) t 


t« l 

t» l 


(18) 


We remark that Eq. ( fl8| ) does not apply to the long time 
scaling of CTRWs, for which it would predict a vanishing 
MSD. In fact, CTRWs do not exhibit a crossover from 


where we define the correlated noise £(i) with (((t)} = 
0 and two-point correlation function: (C(*i)C(^ 2 )) = 
2aK(ti)S(ti — tz). This explicit time dependence clearly 
signals that £(f) is a non stationary noise. It is easily 
shown that the MSD of l*(t) is identical to the one of 
Y(t) given by Eq. (13). However, even if they share the 
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same MSD, Y(t) and Y*(t) provide different PDFs. In¬ 
deed, Y*(t) corresponds to a time rescaled Brownian mo¬ 
tion X(t*) with transformation: 


= / K(r)dT. 

Jo 


( 20 ) 


In the case of the usual Brownian motion the correspond¬ 
ing diffusion equation has a Gaussian solution: P(y, t) = 

for the initial condition P{y, 0) = S(y — 


. (v-v o) 

4crt 


2/o)• Since Y*(t) is just Brownian motion in the rescaled 
time t*, we obtain similarly a Gaussian solution, provided 
we choose the same initial condition: 


P{y,t) = 


1 


\/ Anat* 


(y-y o) 

4 at* 


( 21 ) 


with t* as in Eq. (201. We see that P(y,t) is a solution 


of the diffusion equation: 


^P{y,t) = aK(t)^P{y,t), 


( 22 ) 


with the time dependent diffusion constant: D(t) = 
aK(t). We remark that Eq. ([9| can be recovered from 
these general results by setting 4>(A) = A Q , i.e., K(t) = 
t a ~ 1 /T(a) and t* = t a /T(l + a). However, in order to 
have exact equivalence, we need to neglect the constant 
multiplicative factors in both I\{t) and t* and make the 
following substitution: er —► a a. 


III. LANGEVIN FORMULATION OF 
ANOMALOUS PROCESSES IN PHYSICAL TIME 

A. DEFINITION OF THE NOISE 

We proceed in this section to derive a Langevin descrip¬ 
tion of the process Y(t) defined in Eqs. ( 2apb ) directly 
in physical time. Starting from the explicit integral ex¬ 
pression Eq. (l8|) we can write the following: 


p+oo ps 

Y(t) = V2a S(s — S(t)) / C( r )dT 

= ^i (-| e(i - T(s,) ) \J 0 f(T)dT 

_ p+ oo 

= V2a / 0(i — T(s))£(s) ds, 

Jo 


ds 

(23) 


where the fundamental relation of Eq. (121 is used to ob¬ 


tain the second equality and we then get the third one 
with an integration by parts. We remark that the bound¬ 


ary term [—0(f — T(s)) fg £(t) dr] | Q is zero trivially 
for s = 0, but it vanishes also for s -A +00 because T(s) 
is increasing, thus always being bigger than any fixed 
(and finite) time t. Written as in Eq. (23), Y(t) is a dif¬ 


ferentiable (although in a generalized sense) function of 


time, so that we can take its derivative and obtain the 
equivalent Langevin equation: 


Y(t) = \/2a£(t) 


with the newly defined noise: 


_ p+00 

t(t)= £(s)6(t-T(s))ds, 

Jo 


(24) 


(25) 


whose properties are fully determined by the choice of 
the waiting times’ distribution, i.e., equivalently of the 
function <i>(s) in Eq. If we recall the independence 

of £(s) and r/(s), we can show that £(t) has zero average 
and two point correlation function: 


(f(ii)f(*2)> = K(h)S(t 1 - t 2 ), 


(26) 


with K(t) being specified by Eq. ( |14| . In Laplace space, 
indeed, (?(ti)C(^ 2 )) is an integral of the two point char¬ 
acteristic function of T(s), which can then be computed 

with Eq. ([TIJ: (?(Ai)f(A 2 )^ = $(Al 1 +A2) ■ By making 
its inverse Laplace transform, Eq. (26) follows straight¬ 


forwardly. Consequently, the character of the noise £(i) 
significantly depends on the choice of the function $(s) 


in Eq. (11). Thus, Eq. (24) defines a new Langevin model 
driven by a generalized and typically non Gaussian noise, 
except possibly for particular choices of the memory ker¬ 
nel K(t), which is able to describe free diffusive anoma¬ 
lous processes with arbitrary waiting times distribution 
equivalently to the subordinated Eqs. ( 2apb ). 

We highlight that so far the standard renewal picture 
underlying conventional CTRWs still applies. Eq. (24) is 
essentially the time derivative of Eq. (|8| expressing the 
process in terms of stochastic increments. Differences 
will appear when external forces act on the diffusion pro¬ 
cesses. This case is discussed in Sec. HVl 


B. CHARACTERIZATION OF THE 
MULTIPOINT CORRELATION FUNCTIONS 


The definition in Eq. (25) enables us to derive a 


complete characterization of the multipoint correlation 
structure of £(t). As a preliminary step, we de¬ 
rive the Laplace transform of the multipoint charac¬ 
teristic function of T(s), i.e., Z(t\, si;...; £jv, Sn) = 

fi(tm ~ ViV G N. This is obtained by 

using the definition of Eq. (2b) as: 


N 


z{ Ar, Si; ...; X N ,s N ) = ( ]J e ~ x ™ Jo m \ . ( 27 ) 


Let us first assume an ordering for the sequence of times: 
Si < s 2 <■■■< sn and compute the corresponding 
Eq. (27). If we rearrange the exponent by separating 
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successive intervals, we obtain: 

Z(Xi, si;...; Xn, sn ) = 

„-Aisr //A, v(s' N )ds' N -...-(\ N + ...+\ 1 )/ 0 ” 1 ) 7 (sj)dsj 


= \ 6 
= ^ e -E^[(E 

N-l 


hi A ™)] /* S r + 1 ’n(s'm) ds^ + i ' 

e -( Sm+1 - Sm )<o(E" =r „ +1 A„) ; (28) 


m =0 


where we define So = 0 to simplify the notation and we 
exploited the independence of the increments of T(s) to 
factorize the ensemble average. Furthermore, their sta- 
tionariness together with Eq. (11) is then used to get 


Eq. (28). However, in the general case where no a-priori 


ordering is assumed, we need to consider all the possible 
ordered sequences. We then introduce the group of per¬ 
mutations of N objects Sn, whose elements act on the se¬ 
quence: s = (si,..., Sn)- When we make a permutation 
of s, we obtain a new sequence with permuted indices: 
s' = (s cr (i),..., Scr(Ar)). All the possible orderings of s 
are thus obtained by summing over all the permutations 
in Sn- If we assume that ct(so) = 0, Vct £ Sn, he., the 
initial time is kept fixed by the permutations, and we use 
the result of Eq. (p8|), we derive: 


JV-l 

Z{Xl, Sl, ... , A N, S^) ^ ) | 0 (S(r(m-t-l) X 

(JGSn ttl—0 

x e -(^( m+ i)-M m )) $ (£"= m +i a„ m ) ( 29 ) 

with the ordering of the permuted sequence being speci¬ 
fied by the product of Heaviside functions. By factorizing 
out the first term, we obtain the fundamental result: 


ZiXusn ■ ■ ■; A n ,s n )=Y1 M x (30) 

<5 iv 

x n e(^(m + i)-*»w) e " (8 ' ( ” +rs ’ ( ” ,)$(E ” = ” +1 A ' w) 


m =1 


As an example, we recover the two-point case: 
Z(Xi, si; A 2 , s 2 ) [33]- If we put N = 2 in Eq. (30) 


and we consider the two possible permuted sequences: 
s = (s 1 , s 2 ) and s' = (s 2 , Si), we obtain: 

Z(Ai, * 1 ; A 2 , s 2 ) = B(s 2 - si) e -si*(Ai+A 2 ) e -( S2 - Sl )*(A 2 ) 
+ 0(si - s 2 ) e _S2$ ^ 1+A2 )e _(si_S2)a,(Al) . (31) 
We now use Eq. (|30|) to compute the correlation func¬ 


tions of £(t). Indeed, we obtain from Eq. (251 VN £ N: 

/*+oo 


(£(*l) • • -£(*2JV)) = 


n 


ds„ 


-1 Jo 


.171= 1 


2N 


2N 


n ^) (n - t ( s ™)) > > ( 32 ) 


where the average is factorized due to the independence 
of the noises. If we recall the Wick theorem holding for 
the white noise £(s) [55] : 


/2N \ 

n^) 

\j =1 / 


1 


N 


N2 N 


E IK^(^(2JV-i+i))^(io-(j))) 


(TE.S 2 N 3 =1 
N 


N 2 N 


e n &{3<t( 2N—j+1) t(T (j)) (33) 


creS 2 N j —1 


and we substitute it in Eq. (32), we can derive: 

2 AT /»+oo 


dSn 


<^(*1) - --c(*2 N ))=j^n e ri / 

a£S 2 N Lm—1 J ^ 

N / 2N \ 

x n 6 { s °m-j+i) - md) (n S(ti - T( Si )) 

i=1 \i=i / 

-N /* + oo 


w £ 


N2 N 


(7E.S2N 


n 

.m =1 ’ 


ds, 


cr(m) 


(34) 


AT 


X ^(^ <J ( 2Ar -i+ 1 ) T ( S o-(j)))<5(to-(j) T (scrfj))) j 

with N integrals being solved by using the delta functions 
obtained from (£(si)... £(s2iv)}- ^ we make a Laplace 
transform of Eq. (34), we obtain an expression involving 


Z( Ai, si; A 2 , s 2 ;...; Xn, sn)■ 

= e 

\j=l / 


A/ - 

n 


r +00 


ds„ 


x ^(A 0 .(i)+A 0 .( 2 Ar), si;...; X a ^ N )+X a (N+i), sn) , (35) 


which can thus be further simplified with Eq. ([30]) . By 
substitution and by making a further permutation of the 
indices, we obtain: 


’ 2 N 


n«x))=3vE z 


JV 2 N 


N /»+00 

n / ds CT /( ri 


\j =1 / aeS2NO-'eSN Lm=l * 

AT—1 


xe ^'(1)4 > (E"=1 A m ) [0 (s CT j (m+1) - s CT / (m) ) x (36) 


X g-( s a'( m +l)- s a'(m))' I '(E™= m +l(A ff( ^'(„))+A a(2A r_ CT / ( „ J+1) )) 


where the N integrals can then be solved by making suit¬ 
able changes of variables. This leads to the following 
result for the Laplace transform of even multipoint func¬ 
tions of£(t): 


£(Ai) • ■ .f(A 2 jv)) — 

AT-] 

e n 


1 




AT2^$fri J l 1 A m , CTeS2 


E x ( 3? ) 


ct'g5jv m=l ^ 


(j 2 n=m +1 


cr(<j'(n)) “I” ^a(2AT-o- / (n) + l 


)) 


X 
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We remark that odd multipoint correlation functions are 
zero; indeed, if we make the substitution 2 N —► 2TV +1 in 


Eq. (321, we obtain an expression depending on the odd 
multipoint correlation functions: (£(si)... £(s 2 at+i)), 


which vanish V1V £ N. The corresponding quantities in 
time are derived by making the inverse Laplace trans¬ 
form of Eq. (37), which can be written as a 27V—fold 
convolution: 


JV-l 


(£(tl) ■ ■ ■ £(^2Jv)} = H 8(U +1 - U) *2N t2N-l,t2N) 


JV-l 


g(ti,t2,-..,t2N~i,t2N)=£ 1 < 5Z n 


^<TeS 2 N ct'GSjv m=l ^ (X/n=m+l('^ CT ( CT '( n )) '^o-(2JV-o' , (n) + l))^ 


(38a) 

(38b) 


with K(t) being the memory kernel defined in Eq. (14). 
The set of Eqs. ( |38a||38b 1 can then be used to compute 
all the multipoint correlation functions of £(f) and con¬ 
sequently of Y(t). It is straightforward to recover the 
two point case of Eq. (26), whereas we provide below as 


an example the four point function. First, we need to 
compute Eq. (38b) in time space: 


g(ti,t 2 ,t3,U) = [K(ti)S(t 2 - ti)5(t 3 )5(U) 

+K(ti)S(ti - t 3 )S(t 2 )S{ti) + K(t 2 )S(t 2 - t4,)S{ti)S(t 3 ) 
+K (ti)S(ti — ti)S(t 2 )S(t 3 ) + K(t 2 )5(t 2 — t 3 )S(ti)S(ti) 
+K(t 3 )S(t 3 - f 4 )<5(fi)<5(f 2 )], (39) 


and then solve the ^ 2 ^ 


N=2 


= 6 convolution integrals 

of Eq. (|38aJ) . This can be done explicitly, so that we 
derive: 

(?(*i)?(* 2 )?(* 3 )?(* 4 )) = [iV(min(t 1 ,t 2 ))x 
x K(\ti - t 2 \)5(ti - h)S(t 3 - t 2 ) + A'(min(ti,t 3 ))x 
x K(\ti - t 3 \)6(ti - t 3 )S(t 2 - ti) + /v(min(ti,t 4 ))x 
xK(\ti - U\)S{t 3 - ti)S(U - t 2 )]. (40) 

We verified that the same similar structure of the time 
dependent coefficients is shared by the six point corre¬ 
lation function. Considering the recursive structure ev¬ 
ident from Eqs. ( 38a|38b ), we conjecture the following 
formula for the even correlation functions in time space 
(with t 0 = 0 kept fixed by the permutations): 


' 2 JV 


1 


JV 


n«‘))=jv2»ii n a (j 


cr(2N— m+1) ^<r(m)) ^ 


\j —1 / aeS2Nm=l 

x 53 e(t 0 

ct'eSn 

^ 1))) • 


(m)) t(T(a'(m— 1 ))) ^ 


(41) 


C. COMPARISON WITH THE SCALED BM 

Once the underlying noise structure of the CTRW is 
revealed by Eqs. (38a]|38b||44j, a comparison with the cor- 


r 


responding multipoint correlation functions of the noise 
£(f) of the SBM reveals important common features of 
these two processes. Indeed, the correlation functions of 
£(f) are obtained straightforwardly by using the defini¬ 
tion of Eq. (19) and the Wick theorem in Eq. (33): 


' 2N 


( IT £(*■?) ) _ N 2 N 55 II K ( tcr (m)) x 

\j =1 / ~ o-gS 2 jv m—1 

X S(t a (2N— m+1) ter (m)) ■ (42) 

Odd correlation functions of £(i) are zero as for £(i). As 
an example to better clarify our discussion, we provide 
the four point correlation function: 

(C(tiX(t 2 )C(t 3 )C(U)) = K(t 1 )K(t 2 )6(t 1 — t 3 )S(t 2 -t 4 ) 

+ K(ti)I\(t 3 )5{ti - t 2 )S(t 3 - U) 

+ K{t 2 )K(ti)5{ti - ti)S{t 2 - t 3 ). (43) 

A first remark has to be done when we set N = 2, thus 
studying the two point correlation function. Indeed, this 
is found to be the same for both the noises £(f) and C(t) 
and equal to Eq. ( 26 1, thus explaining why the corre- 


JV 


sponding integrated processes Y(t) and Y*(f) share the 
same MSD. On the contrary, differences are evident only 
if we look at the higher order correlation functions. Thus, 
the two integrated processes are distinguishable only by 
looking at quantities dependent on these higher order cor¬ 
relation functions, e.g., the PDFs or the corresponding 
higher order correlation functions of the integrated pro¬ 
cesses. Furthermore, by comparing Eqs. ( Fipi ) , we can 
observe the same similar structure of the delta functions, 
typical of Gaussian processes, but with a different corre¬ 
lated and mainly not factorizable time structure of the 
coefficients in the case of £(f), which depends on the dif¬ 
ference between successive time in the ordered sequences. 
This ultimately causes its non Gaussian typical charac¬ 
ter. In fact, in the specific case of a constant memory ker¬ 
nel, for all times or in some scaling limit, the two noises 
coincide and reduce to a standard Brownian motion. 



























IV. MODELS WITH EXTERNAL FORCES 


We now consider models of anomalous processes in the 
presence of external forces Basra. Let us first focus 
on the random walk picture of the CTRW and assume 
that external force fields, which depends on the position 
of the walker, only modify its dynamics during the jumps. 
In the continuum limit, these forces are then naturally 
included in the Langevi n eq uation of the process X(s), 
thus modifying Eqs. ( papb| into [31] : 


X(s) = F(X(s)) + V2a£(s), 
T(s) = rj(s), 


(44a) 

(44b) 


where the function F(x) satisfies standard conditions 
1 55] . However, different scenarios may be observed in 
experiments where forces can modify the position of the 
walker also during the waiting times between different 
jumps, without ultimately changing the underlying wait¬ 
ing times distribution. For instance, we would expect 
this situation to occur for the motion of an organelle in¬ 
side the cytoplasm of a cell, which is freely migrating or 
driven by an external field. This different situation turns 
out not to be easily described with the time-change tech¬ 
nique, as it is not clear how to modify the Langevin sub¬ 
ordinated equations in order to take into account these 
further changes in the position variable. However , the 


characterization of the noise £(t) provided by Eqs. (38a 


38b I, or equivalently by Eq. Mil) , enables us to describe 


it by defining a new class of models, defined with the 
Langevin equation: 


Y(t) = F{Y(t)) + V2a£(t). 


(45) 


The difference between the dynamical behaviours gen¬ 
erated by the two models becomes clear when we look at 
their simulated trajectories. In Figure[2]we plot the paths 
of Y(t) obtained both via subordination of Eqs. ( |44a[ 
44b I (panel b) and via integration of Eq. (45) (panel a) 
for a linear viscous-like force F(x) = —'yx with 7 posi¬ 
tive real constant. On the one hand, in the subordinated 
dynamics (dotted arrows, panel b) we observe time inter¬ 
vals where the corresponding anomalous process Y (t) is 
constant, meaning that the walker, in the corresponding 
renewal picture, is waiting for the next jump to occur 
without any force being able to modify its position. On 
the other hand, during these same intervals the process 
Y(t) generated by Eq. (451 is rapidly damped towards 
zero (dotted arrows, panel a), meaning that the walker is 
being driven by the external force. While indeed exter¬ 
nal forces act only during the jump times in the standard 
subordinated case, in our case they affect the dynamics 
of the system for all times, without intrinsically modify¬ 
ing the waiting times distribution and equivalently the 
relation between the number of steps and the physical 
time. We mention that another scenario involving exter¬ 
nal fields directly modifying the waiting time distribution 
of the random walk has been recently discussed in 153 , 



Figure 2. (Colors online) Simulated trajectories of a CTRW 
with a linear viscous-like force acting along its all time evolu¬ 
tion (panel a, Eq. (451) or acting only during the jumps (panel 
b, subordinated Eqs. ( |44a|44b I). Numerical algorithms are 
adapted from [49]. The difference on how the force affects 
the dynamics during trapping events is evident (dotted ar¬ 
rows): (a) the force acts on the particle, thus damping Y(t) 
towards zero; (b) the force does not act, so that the particle 
gets physically stuck and Y(t) is kept constant. 


but this formalism does not have an evident connection 
with ours. 


Clearly, the inclusion of a force changes the renewal 
picture for the position variable Y(t), which can no longer 
be expressed as a superposition of i.i.d. position incre¬ 
ments as in Eqs. These increments now depend on 

the accumulated position up to the time before the jump. 
However, the process T(s), i.e., the stochastic process of 
the jump times parametrized by the arc-length still rep¬ 
resents a renewal process, since the waiting times are 
unaffected by the force. 

In the following, we present a comparison of the MSD 
obtained from Eqs. ( 44a||44b]|45| ) for a tempered stable 
subordinator as in Sec. Ill CT and for different choices of 
the external force F(x). Except when explicitly stated we 
assume zero initial condition, so that the MSD coincides 
with the second order moment. We recall that the model 
of Eq. ( |~i~5| ) defined with the time scaled noise ((t) instead 
of £(f) provides the same MSD. 


A. CONSTANT FORCE CASE 


We first look at the case of a constant homogeneous 
force field: F(Y(t)) = F with F € R + , for which Eq. (45) 
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becomes: 


Y(t) = F+V2a£{t). (46) 

This equation can be solved formally for the exact tra¬ 
jectory of Y(t): 


Y(t) = Ft + \/2a f £(r) dr 

Jo 


(47) 


and then used, together with Eq. (26), to derive the MSD: 

(Y 2 (t)) = F 2 t 2 + 2a f K(t) dr (48) 
Jo 

or equivalently in Laplace transform as a function of <£(s): 


= WT 


2 F 2 2<t 


A 3 A4> (A)' 


(49) 


In the subordinated case, the MSD is computed with the 


same technique of Eq. (13) but with the different variance 
(X 2 (s)) = (F 2 s 2 + 2 as). In Laplace space we obtain: 


F 2 (A)) = 


2 F 2 


+ 


A(3>(A)) 2 A$(A)' 


(50) 


The Laplace inverse transform of both Eqs. p9|[50 1 is 
plotted, together with their corresponding scaling be¬ 
haviours, in Figure[3](main panel and inset respectively). 
In the small time limit, we find that both share the same 


power-law scaling of Eq. (18). However, they differ be¬ 


tween themselves and with Eq. ( |18| when we look at the 
scaling for long times. On the one hand, Eq. (49) pro¬ 
vides the long time scaling: ( Y 2 (t )) ~ F 2 t 2 . Hence, the 
constant force in this limit induces a crossover from subd- 
iffusive to ballistic dynamics. Examples of this nonlinear 
behaviour have been recently discovered in the dynam¬ 
ics of chromosomal loci, which exhibit rapid ballistic ex¬ 
cursions from their fundamental subdiffusive dynamics, 
caused by the viscoelastic properties of the cytoplasm 
[ICES]. Furthermore, it is evident that the exponential 
dumping of the waiting times’ distribution does not affect 
the long time scaling, differently from the corresponding 
scaling of Eqs. (50), which turns out to be (Figure |3j 
inset): 


(r 2 M) 




2F 

T(l+2a) 


t 2 


A* 7^ 0 

fj, = 0 


(51) 


Thus, we find the same crossover to ballistic diffusion 
when /z yA 0, but with different /i-dependent scaling coef¬ 
ficients, whereas in the CTRW case (/z = 0) this crossover 
pattern is lost and the power-law scaling is conserved, al¬ 
though with a different exponent. 



Figure 3. (Colors online) MSD of an anomalous process with 
tempered stable (a = 0.2) distributed waiting times in the 
presence of a constant force acting throughout the all dynam¬ 
ical evolution (main panel) or only during the jump times 
(inset). These two different scenarios are obtained with the 
£—driven process or with the subordination technique, i.e., by 
numerical Laplace inverse transform of Eqs. ( 49|50 1 respec¬ 
tively. The different long time scaling is evident: (main) the 
£—driven process exhibits crossover to ballistic diffusion in all 
cases and without any dependence on the tempering param¬ 
eter /z; (inset) the time-changed process exhibits crossover to 
ballistic diffusion with /i-dependent scaling coefficient when 
/z yf 0, whereas it still scales as a power-law with exponent 
2 a for /.t = 0. 


B. HARMONIC POTENTIAL CASE 


We now consider an external harmonic potential, lead¬ 
ing to a friction-like force: F(Y ( t )) = ~'yY(t) with 7 real 
positive constant. Thus, Eq. (45) provides the following: 


Y(t) = -7 Y(t) + V / 2cr £(f). 


(52) 


As before, we can solve formally Eq. (52) for the trajec¬ 
tory of Y ( t ) and use it together with Eq. (26 1 to compute 
the Laplace transform of the corresponding MSD: 


F 2 (A)) = 


2(7 


(A + 2 7 ) $ (A)' 


(53) 


On the contrary, in the subordinated case we can proceed 
as in Eq. (13) by substituting: (X 2 (s)) = ^ (l — e -27S ). 
One can thus obtain the result below: 


F 2 (A)) = 


A [2 7 + $(A)]' 


(54) 


We plot in Figure |4| the numerical inverse transform of 
Eqs. ( |53p4l ) ( main panel and inset respectively), along 
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with their scaling behaviour for small times. While the 
small time scaling is in both cases the same as in Eq. (18), 


we observe a different behaviour in the long time limit. 
Indeed, we find for Eq. (531 the following scaling laws: 


(Y 2 (t)) 


7 r (°0 


t° 


m ¥= o 

/x = 0 


(55) 


Thus, in the CTRW case the MSD decreases as a power- 
law towards zero. If we recall that this process is equal 
to the SBM up to the MSD, this is the same anomaly al¬ 
ready reported in m- However, we also show that Y(t ) 
correctly converges to a plateau for /i ^ 0, this being 
the expected dynamical behaviour of confined diffusion. 
By recalling that the waiting times are tempered stable 
distributed, the interpretation of the mentioned anomaly 
becomes clear. Indeed, the truncation of the power-law 
tails of the waiting times’ distribution is fundamental to 
let the system find a stationary state, so that the MSD 
can converge to a plateau , which is typical of confined 
diffusion. In fact, no damping of the tails is done in the 
CTRW case, meaning that very long trapping events may 
still happen with non zero, but small probability. Thus, 
if we wait long enough, i.e., in the long time limit, these 
events eventually occur. However, Eq. (52) establishes 


that the system is affected by the external linear force 
also during such events, which then damps all the oscil¬ 
lations of the system. This clearly implies that the MSD 
should decrease to zero, because the system is not able 
to disperse and gets immobilized in Y = 0. On the con¬ 
trary, in the subordinated case the effect of the external 
force is stopped during the trapping events, so that the 
system does not get trapped in the zero position in the 
long time limit. Indeed, the MSD for different values of 
/i share the same long-time plateau: (Y 2 (t )) ~ 


V. CONCLUSION 

In this work, we identified the underlying noise struc¬ 
ture of a free diffusive CTRW with an arbitrary wait¬ 
ing times distribution and we defined its correspond¬ 
ing stochastic force. This enabled us to write a new 
Langevin equation, describing its dynamics directly in 
physical time and equivalently to the original formu¬ 
lation obtained with the subordination technique. We 
then derived a general formula, both in Laplace space 
and in physical time, providing all its multipoint correla¬ 
tion functions, which, although presenting the same time 
structure of Gaussian processes, have time dependent co¬ 
efficients with a non factorizable dependence on the mem¬ 
ory kernel generated by the corresponding subordinator 
of the equivalent time-changed formulation. Thus, ex¬ 
cept for the specific choice of a constant kernel, which 
recovers the factorizability of these coefficients, but re¬ 
duces the noise to a standard Brownian motion, our new 
£—noise was shown to be naturally both non Gaussian 
and non Markov. 



Figure 4. (Colors online) MSD of an anomalous process with 
tempered stable (a = 0.2) distributed waiting times in the 
presence of a linear viscous force acting at all times (main 
panel) or only during the jump times (inset). These two 
cases are obtained with the driven process or with the 
subordination technique, i.e., by numerical Laplace inversion 
of Eqs. ( 53|54 1 respectively. Whereas for small times the 
two processes exhibit subdiffusive scaling, their long time be¬ 
haviour clearly differs: (main) the MSD of the driven pro¬ 
cess decreases to zero in the CTRW case (/x = 0), whereas it 
converges to a /x-dependent plateau for (i ^ 0; (inset) in the 
subordinated case all the curves converge to the same plateau. 


We then investigated the dynamics exhibited by pro¬ 
cesses driven by the £—noise in the presence of external 
force fields and compared it with the one observed for 
usual subordinated processes. In general terms, we found 
that these processes belong to a new class of CTRW-like 
processes where external forces are exerted on the system 
at all times, i.e., both when the corresponding walker 
jumps or waits for the next jump to occur. Clearly, this 
is different from the original subordinated model, where 
external forces are implicitly assumed to modify the dy¬ 
namics only during the jump times. Consequently, dur¬ 
ing the typical trapping events of subdiffusive dynamics 
the anomalous process Y (t) becomes constant on the one 
hand, when it is generated via subordination, or it is de¬ 
terministically driven by the force on the other hand, 
when it is driven by the £—noise in physical time. 

Furthermore, we found that these processes have the 
same MSD of those obtained with the characteristic noise 
of the SBM with time dependent diffusion coefficient 
being a function of their memory kernel. This rela¬ 
tion indeed both provides a better interpretation for the 
anomaly reported in [JO] and show that the correct scal¬ 
ing of the MSD typical of confined motion can be ob- 
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tained by choosing more general kernels, which prevent 
an unbounded decay of the diffusion coefficient. 

For future work it will be interesting to investigate the 
aging and ergodicity breaking properties of our new class 
of processes jSH! ■ This might further differentiate it from 
other anomalous processes such as the SBM. The prop¬ 
erties of time-integrated observables of the £—driven pro¬ 
cesses, which are expressed as functionals of their fluctu¬ 
ating trajectories, are also an open problem. For func¬ 
tionals of CTRWs, closed-form evolution equations can 
be derived that generalize the Feynman-Kac framework 
to anomalous processes ESI EDI- A further generaliza¬ 
tion to anomalous processes with arbitrary waiting time 


distributions has recently been obtained I .'ill , which high¬ 
lights the connection between the waiting time distribu¬ 
tion and the memory kernel appearing in the fractional 
evolution equations. It will be highly interesting to in¬ 
vestigate whether similar closed form equations can be 
formulated for functionals of trajectories driven by our 
£—noise. 
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